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Abstract 

The traditional method of teaching canonical transformations involves the introduction of gener- 
ating functions of various types. This method obscures the underlying structure of the Hamiltonian 
least-action principle, and can make a straightforward concept seem arcane. In this article, I present 
a method for calculating canonical changes of variable in Hamiltonian mechanics using a differential 
approach which is much more straightforward. This method handles canonical variables directly, 
but also returns the correct equations of motion for non-canonical variables. It is also much more 
algebraic than generating functions, making it easier to present in a systematic manner. 
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I. INTRODUCTION 



Canonical transformations are the strongest tool presented by Hamiltonian mechanics 
over Lagrangian mechanics, specifically transformations that yield the action-angle vari- 
ables and carry conserved quantities. There are many equivalent definitions of canonical 
transformations of variables, but they all boil down to the statement that a set of vari- 
ables (P, Q) are canonical transforms of the variables (p, q) if they preserve the symplectic 
structure of the equations of motion. 

The most commonly presented method for obtaining canonical transformation are the 
generating functions outlined in many classical mechanics textbooks (Ofl to name a few). 
Derivatives of these generating functions yield the new set of canonical variables in terms of 
the old set, and in this sense they are straightforward. Most classical mechanics textbooks 
emphasize obtaining the generating function, whence they go about carrying out canonical 
transformations this way. 

In this paper, I present a method of carrying out canonical transformations which is 
equivalent to the generating function method, but which is more flexible in the sense that 
one does not need the generating function to carry out a canonical transformation. This 
has the advantage that new canonical variables can be obtained more or less by algebraic 
methods, and so new canonical variables may be introduced readily to any Hamiltonian 
and any new set of variables, regardless of whether a generating function is known. It 
furthermore has the advantage of generalizing to non-canonical coordinates if one should 
ever need these things. 

It is important to stress that no new physics derives from any choice of coordinate system. 
The use of canonical variables is merely a convenience, although a very useful convenience, 
in describing dynamical systems. 

In Section II, I present a concise derivation of the method that springs from the Legendre 
transformation from Lagrangian to Hamiltonian action integrals. The general rules for 
carrying out canonical transformations using this new method are presented in Section III. 
Section IV presents a series of applications and academic problems which use this method. 
These range from the change of independent variable to a use of complex-valued canonical 
variables to cope with the two dimensional harmonic oscillator potential. Section V presents 
the Dirac "raising" and "lowering" operators from quantum mechanics as an example of non- 



canonical variable transformations in classical mechanics. 



II. LAGRANGIAN TO HAMILTON! AN ACTION INTEGRALS 

The Lagrangian action integral of classical mechanics, with independent variable s, is 
given by 

S = L(q,q,s) ds (l) 

where overdot denotes a derivative with respect to s. The equations of motion are determined 
by the trajectories in configuration space for which the action integral is minimized, viz. 
5S = 0. I will assume the reader is familiar with the calculus of variations enough that this 
first variation statement is clear. If not, there is an excellent general reference by Weinstock 
ij] on the calculus of variations. 

Now, the Hamiltonian is a Legendre transformation on the Lagrangian. The canonical 
momentum is defined by 

9L 
dq l 

for whichever component is desired. The Hamiltonian is then defined by the transformation 

H = Pl q l -L(q,p,s) (3) 
Inserting this new definition into the action integral above gives 

S = J M-H(p,q,s)) ds (4) 
Hamilton's equations are a direct consequence of minimizing this integral with respect to 

(p,q)- 

At this point I am taking liberties with the rigors of differential forms, but for our purposes 
we will not get into too much trouble. Notice that q ds = dq, so now the action integral is 
given by 

S = J p t dq % — H(p, q, s)ds (5) 

where the implied summation is used. This is the action integral of Hamiltonian mechanics. 
The previous version works fine if we want s to be the independent variable, but this differ- 
ential form will be the focus of the remainder of this note. It is instructive to apply the least 
action principle to equation (j3J) to see that this does indeed return the familiar Hamilton 
equations of motion. 



III. CANONICAL TRANSFORMATIONS & GENERATING FUNCTIONS 



The method presented here recasts canonical transformations from a method of gener- 
ating functions to exploiting the differential form of the Hamiltonian action integral. The 
differential form p % dq l — Hds contains all the information necessary to carry out canonical 
(and non-canonical) changes of variable. 

A set of variables (p, q) are canonical if their action integral is in the differential form in 
equation (J5J). Because of the way we constructed the action integral from the Lagrangian, 
we know this statement to be true because we have used the definition of the canonical 
momentum to reach this point. 

Suppose the Lagrangian is expressed in the generalized coordinate sets such that 

L'(Q,Q,s')=L(q,q,s) (6) 

Such a statement clearly models the same dynamical system. The action integral for each 
one is given by 

J L'{Q,Q,s')ds' = J L{q,q,s)ds (7) 
The Hamiltonian action integrals are then given by 



PQ - H'J ds' = J (pq- H) ds (8) 

The differential form is then given by 

J {P dQ- H'ds') = J {pdq- Hds) (9) 

Both integrals lead to Hamiltonian equations of motion, therefore (P, Q) and (p, q) are both 
canonical variables. By construction, these Hamiltonians correspond to the same physical 
system, and therefore the differentials must be equal up to an exact differential which van- 
ishes upon first variation. Therefore, canonical changes of variable preserve the differential 
form 

(P dQ - H'ds') = (pdq- Hds) + (-j-ds) (10) 



\ ds 

where F is the generating function for the canonical transformation just carried out. The 
conventional approach from here would be to place a heavy emphasis on all the ways F 
can be used. This approach has its uses, but is frequently overemphasized to the point 



that many professional physicists rely entirely on generating functions as the characterizing 
feature of canonical transformations. 

An instructor interested in a broader approach should strongly consider looking at empha- 
sizing the full differential relation in equation ( flOj) . and not overemphasize the importance 
of the generating function. The advantage granted here is that one does not need to look 
up a generating function and calculate old variables in terms of new variables and carry out 
an inversion. This utility is perhaps best demonstrated with a series of examples. 



IV. EXAMPLES 



As canonical transformations are a practical tools in classical mechanics, their uses are 
best illustrated by applications. Therefore, it is now prudent to consider a few examples 
of how this method works in practice. The goal of this section is to provide a number 
of examples of increasing sophistication of using this more direct method, bypassing the 
generating function entirely. Several of these examples are well-known and handled with 
generating functions in standard textbooks. The section on harmonic potentials, to the 
knowledge of the author, cannot be found anywhere else. 



A. Changing Independent Variable 

Suppose we wish to make one of the q variables the independent variable. This is simple 
enough. We identify q* = s'. Immediately, we find that H' = — p*, and that (q*)' = s and 
(p*)' = —H by simple inspection of the differential form. 

The benefit here is obvious and intuitive. In the Lie algebraic language of quantum 
mechanics the momentum is the hermitian operator which generates the translation Lie 
group - p x is the generator of x translations. That is made clear here: since x is the new 
independent variable, and we want to translate our other dynamical variables through x, the 
generator of those translations is the new "Hamiltonian" p x . Thus, this formulation keeps 
Hamiltonian mechanics close to its modern formulation of Lie groups (see, e.g., 
an exercise, show that 

df df 



dx dx 



I2|). As 

+ {f,P*} (11) 



where {A, B} is the Poisson bracket. This is equivalent to the above statement of p x gener- 
ating translations in the x-direction. 



B. Constant Velocity Transformation 

In accelerators, the chosen Frenet-Serret coordinates are chosen around a design trajectory 



(see, e.g., js]) for a particle moving at exactly the design energy and exactly where the 
trajectories are intended. In this context, the useful longitudinal coordinates are those 
which deviate from the design particle trajectory. 

Suppose we want to make a new longitudinal variable if) = s — x/v the new independent 
variable, and let x' — x + s/v describe variation from a given longitudinal trajectory. Such 
a coordinate system is of practical importance to many accelerator physics applications. I 
consider only two dimensions here, as the motion transverse to v is left unaffected by this 
transformation to a moving frame. 
Start with 

dip = ds — dx/v = ds' (12a) 
dx' = dx — vds (12b) 
By direct insertion, this gives the differential form being required to satisfy the constraint 

pdx — Hds = p'dx' — H'dip = p'(dx — vds) — H'(ds — dx/v) (13) 

Then we must satisfy the equations 

p = p ' + H'/v (14a) 

-H = -p'v - H' (14b) 
Direct inversion of these equations gives 

P' \ ( P + H/v 



-H' I \ pv-H 



(15) 



Therefore 

H' = H(p', x'; ip) — p'v (16) 

and the equations of motion follow from this transformation, and are symplectic by construc- 
tion. It is also an interesting and straightforward exercise to watch how a time rescaling 
propagates through the action integral. 



C. The Harmonic Oscillator and Action- Angle Variables 



A frequent example of canonical transformations are the action-angle variables, and the 
most common first example of that is the harmonic oscillator. This is usually a student's 
first introduction to generating functions, canonical transformations, and all the formalism 
that follows. It would be nice to use this example as a swift introduction to the utility of 
this differential approach. 

The objective of action-angle variables is to write the Hamiltonian as independent of the 
canonical coordinate. Let us then consider the simple harmonic oscillator 

H = \ (P 2 + c ¥) (17) 

In the differential form, we want to find a set of coordinates which satisfy 

pdq - Hdt = Jd9 - H\J)dt (18) 

where it would be nice (although not generally necessary) to keep the time variable un- 
changed. For a conserved system, given any initial energy E , the particle will trace out the 
ellipse defined by 

E = ±(p 2 +u 2 q 2 ) (19) 

We also know that it traces out these ellipses at a fixed frequency. We may therefore consider 
rewriting the coordinates as 

p = AVJ sin6 (20a) 

x = cos 9 (20b) 

where 9 parameterizes the point on the ellipse and \fj parameterizes the amplitude. Insert- 
ing this into the old Hamiltonian gives 

H'= (A 2 + u 2 B 2 )j (21) 

and the differential takes the form 

A\fl sin# d (b^/J cos 9} = Jd9 (22) 
Expanding the differential and looking at the left-hand side, we find 



AB ( i sin 29 d J + i ( 1 + cos 29) Jd9^j = 
AB (d[J sin29] + ^Jd9 



(23) 



The first term is an exact differential (in fact, it is a representation of the generating function 
for this canonical transformation) and does not contribute to the equations of motion. The 
second term requires that AB = 2 to match with the Jd9 term above. 
Looking at the transformed Hamiltonian, it would be most convenient if 



We have thus obtained the action-angle variables and their corresponding Hamiltonian for 
the harmonic oscillator, and the generating function for the transformation as a byproduct. 

D. Complex Variables and the Two Dimensional Harmonic Oscillator 

As an exercise in this method that is perhaps a bit more interesting, consider a double 
harmonic oscillator with identical frequencies in the x and y directions. This example 
is worked out using generating functions in the underused Corben and Stehle The 
Hamiltonian for this system is given by 



it might be convenient (and useful) to express this in terms of complex variables, since the 
potential may be written as uj^zz* where z = x + iy. This requires that we find the canonical 
momenta conjugate to these complex valued variables. 

It would be very difficult in general to divine the correct generating function for these 
variables. Using the differential approach, however, it is simply a matter of going through 
the process, x = (z + z*)/2 and y = (z — z*)/(2i). Therefore the differential form is given 



A = y/2co^ 



(24a) 




(24b) 



H=\{vl+V 2 y ) + \ul{x 2 + y 2 ) 



(25) 



by 



p x (dz + dz*)/2 — ip y {dz — dz*)/2 — Hds = 

Pi dz + p 2 dz* — H'ds 



(26) 



Pairing off coordinates and again follow through with the algebra we find that 



V = P = iPx ~ iPy)/2 



(27a) 



p 2 =p* = (p x + ip y )/2 



(27b) 



S 



and then, by direct substitution, we get that the new Hamiltonian is given by 

H' = 2pp* + \w\zz* (28) 

which, by a simple rescaling by \/2 (note that in this formulation it is straightforward to see 
that p' = ap and x' = x/a are canonically conjugate) gives the Hamiltonian as 

H' = pp* + ulzz* (29) 

The equations of motion can be immediately determined. As another exercise, determining 
the action-angle variables in this new coordinate system is useful for perturbation theory. 
Try letting 

p = y/iJutoe* (30a) 
z = y/-iJ/uj Q e- ie (30b) 

and their complex conjugates. Using the same process as for the one-dimensional harmonic 
oscillator above, we can see that this leads to the correct differential form. By direct sub- 
stitution into the Hamiltonian, this gives 

H' = VJ7*uJo (31) 

which is a pure function of the actions and is angle-independent. It is interesting to note at 
this point that the relative phase between the x and y oscillations is actually stored in the 
action coordinate and not in the angle coordinate, and furthermore that there is only one 
angle variable for a two-dimensional system. 

This particular set of coordinates is convenient for harmonic potentials, which can be 
written as the real or imaginary part of a power series in the complex plane. It greatly sim- 
plifies canonical perturbation theory for these potentials, as in the example below. Consider, 
for example, the Hamiltonian described above with an additional sextupole term: 

H = \ (Pi + Pi) + \ul (x* + V 2 )+^ (* 3 - 3xy 2 ) (32) 
In our complex coordinate system, this is simply given as 

H' = pp* + ulzz* + ^[z 3 + (z*) 3 ] (33) 
9 



where S' 3 is rescaled as appropriate. In the action-angle variables, we can write this simply 

as 

H' = ^co + [J 3/2 e 139 + (J*) 3/2 e- i3t> ] (34) 

This is algebraically much simpler to work with than the expression in the x — y coordinate 
system, and a good comparison exercise for the student would be to calculate the first order 
correction to the angles using perturbation theory in the x — y variables and the z — z* 
variables and observe which is more convenient. 



E. Harmonically Confined Particle in a Laser Field 



The previous example required the student to work through one canonical change of 
variables. In a typical application, multiple changes of variable may be needed, and the 
current method makes these run by much more smoothly. Problems with fast-oscillating 
and slowly varying coordinates are ubiquitous in plasma and accelerator physics. Examples 



include various ponderomotive forces from ion traps to laser wakefield acceleration 



the free-electron laser instability 12|, |13J], and various forms of electric propulsion 14]. 



Ill, 



In the previous section I considered a harmonic oscillator. The student can go a step fur- 
ther and look at the problem of a harmonically confined particle in a laser field propagating 
in the ^-direction with a vector potential given by 



Ai = A cos [u(y/c — £)] x 



(35) 



The Hamiltonian, which is effectively two dimensional, is then given by 

H = \ { ( Px ~ ~c A ° cos [w(z//c ~^) 2 +p y} + 



(36) 



-u 2 (x 2 + y 2 ) 



This can be solved fairly directly, using multiple canonical transformations. First, expand 
the Hamiltonian out and write it in terms of the usual x — y action-angle variables derived 
above to obtain 



H = (J X + Jy)U!Q — 



(a/ J x UJq COS 9 x j ^AqCOS 



oj I J J y j uq sin 8 y — t 



(37) 
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where I have assumed that eA /c <^ y/ J x uq and that higher order term can be ignored. We 
can rearrange the perturbing term as 



cos 



sm 



u 
c 
u 
c 



cos 9 X cos 
\J Jy/ouo sin y 
yj Jy/^o sin 6 y 



^ ^yj Jy/ujo sin y - tj 

COs(9 x — Cdt) + COs(^ + Ut)) + 

(sm(9 x + ut) — sm(9 x — ut)) 



(38) 



Our first new change of canonical variables is to select the slow oscillations by choosing 
'&x = 9 X — ut. Then terms that go as e^+ 2 ^ are rapidly oscillating and can be dropped on 
the average. This is commonly referred to as a transformation to the rotating coordinate 
system. 

To accomplish this canonically, look at the action integral, which in the action-angle 
variables looks like 

S = J J x d9 x + J y d9 y - H'dt (39) 
since we are moving into the slow-varying coordinates, we want 



J x d9 x + Jyd9 y - H'dt = J x dd x + J y d9 y - Udt 



(40) 



and need to solve for H and J x . J y and 9 y remain unchanged. In the differential form, this 
criterion reads: 

J x d9 x - H'dt = J x d9 x - Udt - J x udt (41) 
and so J x = J x and T-L — H' — J x u. This new, new Hamiltonian is given by 



le 
2c 



y/ J x u A 1 cos 



V. = J X (0J ~0j) + JyU + 

cos-# T — 



sm 



^yj 'jy/u sin 9y 
-\Jj y /u sin 9y 



(42) 



sin i? 



This is useful because we have used canonical transformations to make the Hamiltonian 
approximately time-independent in the slow-varying frame. So long as w > w this Hamil- 
tonian is conserved, which is a large improvement to the equations of motion. Treating A as 
an expansion parameter, it is then possible to obtain a perturbation expansion in these new 
canonical variables, which were obtained directly with no intermediate generating function. 
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V. NON-CANONICAL VARIABLES 



Canonical variables are nice because the dynamics are known to be symplectic, which 
encodes the conservation laws into the dynamics. However, there may come a point where 
non-canonical variables may be useful. The generating function by design cannot cope 
with this. However, this method will still yield the correct equations of motion, even if the 
equations of motion are not symplectic. This method has been used to study magnetohydro- 



dynamics [15H17] where the use of Eulerian variables is preferable, even though they are not 
canonical. Presented here is an example of such variables familiar to most undergraduates. 

For the two-dimensional harmonic oscillator, we introduced z = x + xy to reduce the 
complexity of the Hamiltonian. We may do something similar for the one-dimensional case. 
As an example, consider the one-dimensional harmonic oscillator 

H = \ [f + c! 2 ) (43) 

Define the non-canonical variables 

a = (p + iq)/V2 (44a) 

a* = {p-tq)/V2 (44b) 

This change of variables might be familiar to anyone who took a quantum mechanics course. 
The Hamiltonian is clearly written as H = aa*. But what are the equations of motion for 
a and a*? We can turn to the action integral and minimize it using the Euler-Lagrange 
equations to obtain the non-canonical equations of motion. We know that p = (a + a*)/-\/2 
and q = —i(a — a*)/v / 2 and therefore the action integral is given by 

% % 

--(a + a*)a-\ — (a + a*)a* — aa* ds (45) 

This does not have the same structure of pdq — Hdt that the canonical transformations 
preserve. Specifically, the sign difference between the coefficient of da and da* prevents the 
differential above from being written in the canonical form to within a total derivative. 

Because the change of variables occurs within the Hamiltonian action integral, we can 
still obtain equations of motion. Minimizing this action integral using the Euler-Lagrange 
equation 

e- f a , d OF dF 

5 / F(x, x')ds = — — - — = 
J ds ox' ox 

12 



gives the equations of motion for the non-canonical variables as 

a* = la (46a) 

a = -la* (46b) 

This is identical to the standard equations of motion one might obtain in the Heisenberg 
picture using these as operators for the quantum mechanical harmonic oscillator. By direct 
addition of these equations we duplicate the standard equation p = —q, then by subtraction 
we obtain that q = p. 

VI. CONCLUSION 

It is the purpose of this paper to present a short series of examples which exploit the differ- 
ential form of Hamilton's action integral to find canonical transformations. The generating 
function is frequently of little direct use aside from finding the new canonical variables, and 
it is usually the emphasis of conventional classical mechanics textbooks to approach canon- 
ical transformations in terms of generating functions. This frequently adds an additional, 
more complicated step to any calculation that is redundant in almost every case. 

By providing a direct derivation of Hamilton's action integral, then working through a 
series of examples that illustrate the usefulness of this method, it is my hope that this 
approach to canonical transformations might take root for its practicality and intuit iveness. 
At the very least, it makes clear where generating functions might come from and how 
to derive them, although from the point of view described here they are not even strictly 
necessary. 
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